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Abundances of Fe, Si, Ni, Ti, Na, Mg, Cu, Zn, Mn, Cr and Ca in the atmosphere of the K-dwarf HD 77338 are 
determined and discussed. HD 77338 hosts a hot Uranus-like planet and is currently the most metal-rich single star 
to host any planet. Determination of abundances was carried out in the framework of a self-consistent approach 
developed by Pavlenko et al. (2012). Abundances were computed iteratively by the program ABEL8, and the 
process converged after 4 iterations. We find that most elements follow the iron abundance, however some of the 
iron peak elements are found to be over-abundant in this star. 
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1 Introduction 

Determining the chemical composition of stars is one of the primary goals of astrophysics. Such inves¬ 
tigations help us to better understand the chemical enrichment of the Galaxy and to make some assumptions 
about the mechanisms involved in element evolution in the interstellar medium, and in stellar atmospheres 
in particular |12| . While studying the Sun, the problem of the abundances of certain atoms necessitated 
a model to explain this. It was finally explained with the introduction of the pp- and CNO- cycles in the 
interior of the Sun. But this was not enough to explain the presence of large amounts of helium. The next 
step in studying the evolution of elements was the introduction of nucleosynthesis theory. Modern scientific 
understanding is that chemical elements were formed as a result of the processes occurring in stars, leading to 
evolutionary changes of their physical conditions. Therefore, the problem of nuclide formation is also closely 
related to the issue of the evolution of stars and planetary system. Recently Jenkins et al. |5] announced 
the discovery of a low-mass planet orbiting the super HD77338 as part of our ongoing Calan-Hertfordshire 
Extrasolar Planet Search [6j. The best-fit planet solution has an orbital period of 5.7361 ± 0.0015 days and 
with a radial velocity semi-amplitude of only 5.96±1.74 ms“^, giving a minimum mass of ('(g'g M^. The 
best-fit eccentricity from this solution is O.OOlgQb and is in the agreement with results of a Bayesian analysis 
and a periodogram analysis. 

According to modern theory, the formation of the nucleus of chemical elements from carbon to iron is 
the result of thermonuclear reactions involving He, C, O, Ne and Si in stars. After the depletion of hydrogen 
reserves, a star’s core starts running a 3 q; reaction, where it produces a number of elements as a result of the 
following transformations: 3'^He —> -|- ^He —-|- 7, -|- ^He —^°Ne -|- 7. 

After reaching a specific threshold temperature, carbon begins fusing with the formation of Ne, Na 


and Mg: 

12c -b 12 c — 

A 20Ne -b 'iHe -b 4, 62MeV, 


12c -b 12c — 

A 23Na -b p -b 2, 24MeV, 


12c -b 12c — 

A 24]y[g ^ _ 2,60MeV. 


Aluminum 

can then be produced by: 24Mg -bp — 

-A 25Al + 7 


The combustion reaction of oxygen is a dual-channel process and causes the presence of Al, S, P, Si 
and Mg. One of the channels is: 

I60 + I60 ^ 30gi + iH + iH + 0,39MeV, 

I60 + I60 —^ 24]y[g 4jjg 4jje _ 0,39MeV, 
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Table 1: Simbad’s list of previous assessments of abundances in the atmosphere of HD 77338 


Teff 

logg 

[Fe/H] 

CompStar 

Reference 

5300 

4.30 

0.36 

Sun 

m 

5290 

4.90 

0.22 

Sun 

m 

5290 

4.60 

0.30 

Sun 

nni 


160 + ^ 27^1 + 4jje + iH - l,99MeV, 

With continuous temperature growth, silicon burning is initiated. This process is described by a 
number of reactions. As a result we can receive, for example, Ar, Ni, S, etc. ®®Ni, after two /3 decays, turns 
into ^®Fe. It is the final stage of the fusion of nuclides in massive stars, which forms the nucleus of the iron 
group. 

The production of heavy elements is provided by other mechanisms. They are called s- and r- processes. 

• s-process or slow neutron capture: formation of heavier nuclei by lighter nuclei through successive 
neutron capture. The original element in the s-process is ^®Fe. The reaction chain ends with It 

is thought that s-processes occur mostly in stars on the asymptotic giant branch. For the s-process to 
run, an important condition is the ability to produce neutrons. The main neutron source reactions are: 

-|- '^He —)• -|- n, -|- ^He —25]y[g _|_ Example of s-process reactions are: 

Elements heavier than H and He are usually called metals in astrophysics. Their concentration is signifi¬ 
cantly less, relative to hydrogen and helium, but they are the source of thousands of spectral lines originating 
from a star’s atmosphere. The abundance of iron depends on a stars age and on its position in the galaxy 
[9]. Metal-rich stars are also known to be rich in orbiting giant exoplanets. High metallicity appears to be a 
major ingredient in the formation of planets through core accretion j5|. 

HD 77338 is one of the most metal-rich stars in the sample of |3] and in the local Solar neighborhood in 
general. Its spectral type is given as KOIV in the Hipparcos Main Catalogue |l7j . However, HD 77338 is not 
a subgiant, as labeled in Hipparcos [5], its mass and radius are smaller than the Sun’s: M = 0.93 ± 0.05 Mq, 
R = 0.88 ± 0.04 Rq. a parallax of 24.54 ± 1.06 mas for HD 77338 means the star is located at a distance 
of 40.75 ± 1.76 pc. Its effective temperature and surface gravity were found T^ff — 5370 ± 80 K, log^ = 
4.52 ± 0.06 [5]. More stellar parameters for HD 77338 and detailed information about its planetary system 
are in [5]. 

Using the Simbad database one can find information on the previous assessments of abundances in the 
atmosphere of HD 77338 (see Tabled]). In most cases the authors only determine the metallicity of the star, 
i.e. the iron abundance. In turn, we recomputed the abundances of many elements which show significant 
absorption lines in the observed spectrum of HD 77338 . 


2 The observations 

The observations of HD 77338 were carried out as part of the Calan-Hertfordshire Extrasolar Planet Search 
(CHEPS) program jTj. The main aim of the program is monitoring a sample of metal-rich stars in the 
southern hemisphere to search for short period planets that have a high probability to transit their host 
stars, along with improving the existing statistics for planets orbiting solar-type and metal-rich stars. The 
high-S/N (> 50) and high-resolution {R = 100000) spectrum of HD 77338 , observed with the HARPS 
spectrograph [10] . was reduced using the standard automated HARPS pipeline and analyzed in this work in 
order to determine the chemical abundances and other physical parameters of the stellar atmosphere. 


3 The procedure 

Firstly, we selected “good” absorption lines for all elements of interest that are present in spectra of the Sun 
and HD 77338. These lines should be not be blended (see 0) and be intense enough in both spectra. We 
selected lists of lines of each element that were to be used for the abundance investigations. We used line list 
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data, which was taken from database of atomic absorption spectra VALD |3 , to compute synthetic spectra of 
the Sun for a plane-parallel model atmosphere with parameters Teg/logg/fFe /H] = 5777/4.44/0.0 |15) . The 
model atmosphere was used to compute the synthetic spectra using WITA6 |13) . building a grid of models 
with different microturbulent velocities 4/= 0-3 km/s with a step size of 0.25 km/s. The shape of the line 
absorption profiles were constructed as Voigt function profiles H(a,v), and a classical approach was used to 
compute the damping effects |20| . To compute the rotational prohle we followed the procedure described in 
Gray [2|. 

All abundance determinations were performed by the ABELS program |14| . Details of the full procedure 
we used is described in |16| . see |3] also for more details on the line selection and fitting procedure. 


4 Results 

4.1 The Sun 

The solar spectrum is well-studied and abundances for the Sun are known to very high accuracy, therefore 
it represents a very good template. Fig. [1] and Fig. [2] illustrate the presence of spectral lines of Cr in the 
observed spectrum of the Sun as a star [8j. Arrows on the plot show the spectral range which was selected 
to compute profiles of two Cr I lines to be used later by ABELS m in the determination of the abundance 
of chromium. We employ a similar selection in the solar spectrum Sun and the spectrum of HD 77338 for 
lines of Fe I, Si I, Ni I, Ti I, Na I, Mg I, Cu I, Zn I, Mn I, Cr I, A1 I, and Ca 1. 

We verified whether our input data are good enough to reproduce the abundances in the atmosphere 
of the Sun. We computed abundances for the Sun using the fits of the theoretical spectra to the profiles 
of the selected lines. In that way we can test our method and estimate the accuracy of our abundance 
determination. Then we investigated the dependence Ea = da/dE", where a and E" are the iron abundance 
and excitation potential of the correspondent radiative transition forming the absorption line. Best fits of 
the selected lines of Fe I in the computed spectra, when compared to their observed profiles in the solar 
spectrum, provides the min Ea of 47=0.75 km/s. The abundances of iron and other elements were then 
obtained using this adopted value for the microturbulence, the results are shown in Table [2j It is worth 
noting that our abundances agree with the reference values within an accuracy of ±0.1 dex. 

4.2 HD 77338 

The model atmosphere for HD 77338 was computed using the parameters determined by Jenkins et al. [5] 
using the SAM12 program |15| . Again, as the first step of our analysis we determined the microturbulent 
velocity in the atmosphere of HD 77338 . The minimum of the slope of Ea provides 4/= 0.75 km/s for log 
N(Fe) = — 4.120 ± 0.07 or [Fe/H] = 0.281 (iteration 1). For other elements we used the same value (4/ = 
0.75 km/s). 

We carried out 4 iterations to determine all abundances. In each next step the ab undances from the 
previous determination were used to recompute the model atmosphere by SAM12 m and the synthetic 
spectra. Each time we are approaching self-consistency by computing the model atmosphere that relates to 
the final metallicity of the star. 

In Table [2] we present our results for 12 different ionic species. We compare our abundances with the 
solar values, obtained using a model atmosphere of 5777/4.44/0.00. They are in good agreement with each 
other. Fig. [2] and Fig. [3] show the line profiles of Cr and Mg calculated using a 4/ = 0.75 km/s. 

In Fig. |4]we show the dependence of [X/H] on atomic number of each element for every iteration. The 
presence of errors can be explained by the presence of noise in the selected spectral lines, along with only 
having a small number of lines to work with for some elements. In Fig. [5] we present the dependence of 
[X/Fe] for the final iteration, where different elements are shown using different plotting shapes depending 
on their mechanism of formation. 


5 Discussion 

We determined abundances for 12 ionic species in the atmosphere of the Sun and the metal-rich exoplanet 
host star HD 77338. Our values for the solar abundances are in good agreement with results from previous 
authors, proving the validity of our method. We used the solar spectrum as a reference to select the proper 
list of absorption lines to be used later in the analysis of the HD 77338 spectrum. 
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Table 2: Abundances in the atmosphere of HD 77338, iteration 4 


Iron 

log N(X) 

log N(X)q ABELS 

log N(X)q 

[X/H] 

[X/Fe] 

usini (km/s) 

Ni 

Al I 

- 5.403 

± 

0.000 

- 5.767 

± 

0.033 

- 5.551 

-b 

0.148 

— 

0.095 

2.33 ± 0.44 

3 

Ca I 

- 5.376 

± 

0.029 

- 5.588 

± 

0.023 

- 5.661 

-b 

0.285 

-b 

0.042 

1.25 ± 0.13 

14 

Cr I 

- 6.085 

± 

0.032 

- 6.345 

± 

0.026 

- 6.441 

-b 

0.356 

-b 

0.113 

1.91 ± 0.14 

22 

Cu I 

- 7.670 

± 

0.058 

- 8.133 

± 

0.067 

- 7.941 

-b 

0.271 

— 

0.028 

2.17 ± 0.44 

3 

Fe I 

- 4.158 

± 

0.038 

- 4.439 

± 

0.023 

- 4.401 

-b 

0.243 

-b 

0.000 

2.06 ± 0.11 

27 

Mg I 

- 4.228 

± 

0.058 

- 4.367 

± 

0.095 

- 4.441 

-b 

0.213 

— 

0.030 

1.50 ± 0.50 

3 

Mn I 

- 5.957 

± 

0.097 

- 6.600 

± 

0.046 

- 6.641 

-b 

0.684 

-b 

0.441 

2.69 ± 0.21 

8 

Na I 

- 5.387 

± 

0.048 

- 5.789 

± 

0.054 

- 5.721 

-b 

0.334 

-b 

0.091 

1.94 ± 0.31 

9 

Ni I 

- 5.367 

± 

0.033 

- 5.756 

± 

0.027 

- 5.821 

-b 

0.454 

-b 

0.211 

2.29 ± 0.15 

17 

Si I 

- 4.111 

± 

0.054 

- 4.469 

± 

0.058 

- 4.401 

-b 

0.290 

-b 

0.047 

2.39 ± 0.13 

23 

Ti I 

- 6.897 

± 

0.040 

- 7.064 

± 

0.028 

- 6.981 

-b 

0.084 

— 

0.159 

1.88 ± 0.13 

24 

Zn I 

- 7.028 

± 

0.065 

- 7.375 

± 

0.048 

- 7.441 

-b 

0.413 

-b 

0.170 

2.25 ± 0.25 
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Our [X/H] correlates well with the condensation temperature of the ions (Tcond), see discussion in |11) . 
This may indicate the presence of a common shell (in the past) and can be an additional criterion for the 
existence of a planetary system around metal-rich stars. 

We also computed usini for both stars. It is worth noting that we determined all parameters in the 
framework of a fully self-consistent approach (see |15| for more details). In general, lines of Mn, Cu can 
not be used to obtain usini because these lines usually have several close components, but in our case the 
parameter v sini was used to adjust theoretical profiles to get the proper fits to the observed special features. 
We believe that fits to Fe I lines provide reasonable measures of the rotational velocity. 

Our results show that the abundances of most elements in the atmosphere can be described well by the 
overall metallicity. However, we found an overabundance of some of the iron peak elements (e.g. Mn, Cu). 
Interestingly, Cu is an element formed through the s-process and its abundance follows that of Fe, whereas Zn 
and elements formed through the p-process, e.g. Ni (see http://www.mao.kiev.ua/staff/yp/TXT/prs.png), 
show a noticeable overabundance compared to iron. It would be interesting to compare these results for HD 
77338 with other metal rich stars to see if this is a common trend for super metal-rich stars. We plan to 
investigate this issue in a following paper (Ivanyuk et al. 2013, in preparation). 
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Fig. 1: The dotted line represents the observed spectrum of the Sun, the solid line is the 
observed spectrum of HD 77338, the dashed line shows the synthetic spectrum computed 
by Wita618 with log N(Cr)= — 6.11 for a model atmosphere of 5315/4.39. This plot was 
used to detect "clean" parts of Cr line prohles marked here by arrows 



Fig. 2: Cr I absorption line profiles computed for a model atmosphere of 5315/4.39 and 
log N(Cr) = — 6.09 using ABEL8 and htted to the observed spectrum of HD 77338 
shown by dashed and solid lines, respectively. The vertical lines show the htted parts of 
Cr I line prohles. 
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Wavelength (A) 


Fig. 3: Computed by ABELS and observed line pro¬ 
files of Mg in the spectrum of HD 77338 with a model 
atmosphere of 5315/4.39, log N(Mg) = — 4.23 shown 
by dashed and solid lines, respectively. The vertical 
lines show the fitted part of Mg I line profile. 



Atomic Number 

Fig. 4: Dependence of [X/H] on atomic number of each element for HD 77338. Open 
circles show values found in the first iteration, filled circles are for iteration 2, open 
squares are for iteration 3, and hnally the stars show the results for iteration 4. 
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Atomic Number 

Fig. 5: Dependence of [X/Fe] on the atomic number of each element. The different 
plotting shapes represent the different formation mechanism of each element. Open 
circles are for a-elements, diamonds show the thermonuclear elemental production, and 
s-process is shown by triangles. 




